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Q ; Abstract 



We study lepton number violating (LNV) process of (/x , e^) conversion in nuclei mediated by the 



. exchange of light and heavy Majorana neutrinos. Nuclear structure calculations have been carried 

^r] ! out for the case of experimentally interesting nucleus ^^Ti in the framework of renormalized proton- 

o 

Q^ ' neutron Quasiparticle Random Phase Approximation. We demonstrate that the imaginary part of 

|: .„.,....„-..„,„™. 

the (/i , e^) conversion rate. This specific feature is absent in the allied case of Oz//3/3 decay. Using 
the present neutrino oscillations, tritium beta decay, accelerator and cosmological data we derived 
the limits on the effective masses of light {m)^^ and heavy (M^^)^e neutrinos. The expected 
rates of nuclear e"*") conversion, corresponding to these limits, were found to be so small that 
^ I even within a distant future the e"*") conversion experiments will hardly be able to detect the 

neutrino signal. Therefore, searches for this LNV process can only rely on the presence of certain 
physics beyond the trivial extension of the Standard Model by inclusion of massive Majorana 
neutrinos. 
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I. INTRODUCTION 



Lepton number (L) conservation is one of the most obscure sides of the Standard Model 
(SM) not supported by an underlying principle and following from an accidental interplay 
between gauge symmetry and field content. Any deviation from the SM structure may 
introduce L non-conservation (LNV). Over the years the possibility of lepton number non- 
conservation has been attracting a great deal of theoretical and experimental efforts since 
any positive experimental signal of LNV would point to physics beyond the SM. The simplest 
extension of the SM allowing LNV processes implies inclusion of massive Majorana neutrinos 



however, the role 
ill consists in the 



with the AL = 2 mass term introducing the necessary source of LNV. 
of neutrinos in LNV processes is more intricate. The fundamental fact 
following: observation of any LNV process would prove that neutrinos are massive Majorana 
particles. This is true even if their direct contribution to this process is negligible and the 
dominant contribution has nothing to do with neutrinos. 

Recent neutrino oscillation experiments established the presence of small non-zero neu- 
trino masses, the fact that itself points to physics beyond the SM. However neutrino oscil- 
lations are not sensitive to the nature of neutrinos: they could be either Majorana or Dirac 
particles leading to the same oscillation observables. 

The principal question if neutrinos are Majorana or Dirac particles can be answered only 
by searching for LNV processes which, as commented above, are intimately related to the 
nature of neutrinos. 

Various LNV processes have been discussed in the literature in this respect (for review 
see j2|). In principle, they can probe Majorana neutrino contribution and provide infor- 
mation on the so called effective masses {mu)a/3 and {M^^)af3 of light and heavy Majorana 
neutrinos (for definition see Sect. |n|). These quantities under certain assumptions are related 
to the entries of the Majorana neutrino mass matrix . 

Among these processes there are a few LNV nuclear processes having prospects for ex- 
perimental searches: neutrinoless double beta decay (OuPP), muon to positron (/i",e"'") 
conversion and, probably, muon to antimuon {fi~,fi~^) conversion j^Q]. 

Currently the most sensitive experiments intended to distinguish the Ma iorana nature o 
neutrinos are those searching for Oz//5/3-decay j^O,!?,!^- The nuclear theory side 
of this process has been significantly improved in the last decade (see also 
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references therein) allowing reliable extraction of fundamental particle physics parameters 
from experimental data. 

The e"*") conversion is another LNV nuclear process searched for experimentally. The 
important role of muon as a test particle for new physics beyond the SM has been recognized 
long time ago. When negative muons penetrate into matter they can be trapped to atomic 
orbits. Then the bound muon may disappear either decaying into one electron and two 
neutrinos or being captured by the nucleus, i.e., due to ordinary muon capture. These two 
processes, conserving both total lepton number and lepton flavors, are the SM processes and 
have been well studied both theoretically and experimentally. The physics beyond the SM 
resides in yet non-observed channels of muon capture: muon-electron a~^e~) and muon- 

n n n r — — — — ^ „ „ „ 

asitron {fi , e"*") conversions in nuclei [15|, [la llTl . 
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(1) 



iA,Z)+fx^-^e- + iA,Zy, 
{A,Z)+fi^^e+ + {A,Z-2y. 

Apparently, the (/x~,e~) conversion process violate lepton flavor Lf and conserve the total 
lepton number L, while {fi^,e^) conversions violate both of them. Additional differences 
between the (/i~,e~) and (/i~,e"'") lie on the nuclear physics side. The first process can 
proceed on one nucleon of the participating nucleus while the second process involves two 
nucleons as dictated by charge conservation 0|- Note also that the (/i^,e~) conversion 
amplitude is quadratic and {fi~,e~^) amplitude linear in the light neutrino mass. Thus the 
second process looks more sensitive to the light neutrino masses. 

The currently best experimental limit on (/i~,e'^) conversion branching ratio has been 
established at PSI 31] for the ^^Ti nuclear target 

r(/i-+48Ti^e+ + ^8Ca 



< 4.3 X 10"^^ 



(2) 



r(/i- + 48Ti ^ Z/^ + 48SC 

Now it is expected a significant improvement of this limit in the near future experiments: 
SINDRUM II (PSI) with ^^Ti target 3ll, MECO (Brookhaven) with ^^Al target [32] and 
PRIME (Tokyo) with ^^Ti target 



In the present paper we study light and heavy Majorana neutrino exchange mechanisms 
of the (/i~,e"'") conversion which are conceptually most natural and simple. One of the 
main motivations of this study comes from the nuclear physics side of this process: the 
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nuclear theory of e"*") conversion is not yet well elaborated and may show new interesting 
features absent in the other LNV processes such as the Oz//3/3-decay. For instance, as we will 
demonstrate, the imaginary part of the (/i~,e^) conversion amplitude in the case of light 
Majorana exchange gives an appreciable contribution to the rate of this process, the fact 
which has not been recognized for a long time. Studying the most simple case of (/x^,e^) 
conversion via Majorana neutrino exchange, we have in mind that this process may receive 
contribution from other mechanisms offered by various models beyond the SM such as the 
R-parity violating supersymmetric models, the leptoquark extensions of the SM etc. Some of 
these mechanisms may involve light or heavy neutrino exchange and, therefore, in the part of 
nuclear structure calculations they may resemble the ordinary neutrino mechanisms. Thus 
our present study can be viewed as a step towards a more general description of (/i~,e"'") 
conversion including all the possible mechanisms. 

Below, we develop a detailed nuclear structure theory for the light and heavy neutrino 
exchange mechanisms of this process on the basis of the nuclear proton- neutron renormal- 
ized Quasiparticle Random Phase Approximation (pn-QRPA) . We calculate the 

nuclear matrix elements of {u_, e^) conversion in ^'^Ti, which serves as target nucleus in 
SINDRUM [31] and PRIME ^ experiments. 

Existing limits on neutrino masses and mixing from neutrino oscillation phenomenology 
and other observational data allow us to estimate typical rate of this process, assuming 
the dominance of light or heavy Majorana neutrino exchange mechanisms. Extremely low 
values for these rates, derived in this way, leave no chance to detect a neutrino signal in 
the {fi~,e^) conversion even within a distant future and, thus, to derive information on 
the effective masses {my)^^ and (M^^)^e from this process. This conclusion, nevertheless, 
does not diminish the importance of experiments searching for (/i~,e"'") conversion since its 
observation would be unambiguous signal of a non-trivial physics beyond the SM. 

The paper is organized as follows. In Sect. |n]we discuss some general issues of Majo- 
rana neutrinos for LNV processes. Sect. IIIII deals with the current limits on the effective 
Majorana neutrino masses entering to the (/i~,e"'") conversion amplitude. The amplitude 
and rate of (/i^,e+) conversion are derived in Sect. IIVI The details of nuclear calculations 
for (/i^,e+) conversion in ^'^Ti are given in Sect. In Sect. IVII we discuss the possible 
impact of (yU~,e"*") conversion experiments on neutrino physics and visa versa. In Sect. IVIll 
we summarize our results and conclusions. 
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II. MAJORANA NEUTRINOS IN LNV PROCESSES 



The finite masses of neutrinos are tightly related to the problem of lepton flavor / number 
violation. The Dirac, Majorana and Dirac-Majorana neutrino mass terms in the Lagrangian 
offer different neutrino mixing schemes and allow various lepton number/flavor violating 
processes jsG, 37, Ssl ]. 

Let us consider the generic case of neutrino field contents with the three left-handed 
weak doublet neutrinos z/^j = (z/^g, z/^^, z/^^) and n species of the SM singlet right-handed 
neutrinos v'^^ = (z^^i, ...z/^„). The mass term for this set of fields can be written in a general 
form as 



3+n 

rrLiU^iUi + H.c. (3) 



^ 3+n 
i=l 

Here Aii, -M-r are 3x3 and n x n symmetric Majorana mass matrices, M^, is 3 x n Dirac type 
matrix. Rotating the neutrino mass matrix by the unitary transformation to the diagonal 
form 

U^M^'^^U = Diag{mi} (4) 

we end up with n + 3 Majorana neutrinos z/j = U^^u'f. with the masses rrii. In special cases 
there may appear among them pairs with masses degenerate in absolute values. Each of these 
pairs can be collected into a Dirac neutrino field. This situation corresponds to conservation 
of certain lepton numbers assigned to these Dirac fields. 

The considered generic model must contain at least three observable light neutrinos while 
the other states may be of arbitrary mass. In particular, they may include intermediate and 
heavy mass states. Presence or absence of these neutrino states is a question for experimental 
searches. 

The favored neutrino model has to accommodate modern neutrino phenomenology in 
a natural way, in particular, to answer the question of the smallness of neutrino masses 
compared to the charged lepton ones. The most prominent guiding principle in this problem 
is the see-saw mechanism. It suggests that the typical scale of matrix elements in Eq. 
(121) is comparable with the masses of charged leptons meanwhile the J^r is associated to 



a large hypothetical scale of lepton number violation like M/^jvy ^ 10^^ GeV. Then the 
diagonalization in Eq. Q brings very light z/^ and very heavy A^^ Majorana neutrinos. 
This mechanism can be realized in various models beyond the SM with significantly lower 
scales, MiNV ~ 1 TeV, leading to the neutrino masses and mixing consistent with the 
observational data. A particular example is given by the class of supersymmetric model 
with bilinear R-parity violation (see, for instance, Ref. js^ and references therein). In these 
models the heavy Majorana neutrinos have moderately large masses ~ ITeV and even lower 
giving them phenomenological significance via a priori non-negligible contributions to LNV 
processes. In the present paper we examine the contributions of light and heavy Majorana 
neutrinos to (/i^,e^) conversion. 

In general, the fiavor neutrino states are the superpositions of light (z/jt) and heavy {N^) 
Majorana mass eigenstates: 



with the masses and respectively. Here U is neutrino mixing matrix. 

Now let us consider LNV processes with two charged (anti-)leptons {la)la, {h)h 



Assume that the characteristic energy scale of this process is go and that light and heavy 
neutrino masses satisfy the conditions: 



Then neutrino contribution to its amplitude Aap can be represented in the form (for more 
details see, for instance, Ref. [4^ ) 




(5) 



initial/final state or with one (/q)Zq, in the initial and another Ip, {Ip) in the final state. 



mk <^ go for V/c, and ^ go for Vfc. 



(6) 




(7) 



where G^, Gn are the corresponding structure factors and 




(8) 



fc=light 




(9) 



are the effective light and heavy neutrino masses respectively. 
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The following comment is in order. If the mixing of heavy neutrino states to the active 
flavors is negligible, the light neutrino sector can be characterized by the effective light 
neutrino mass matrix Ai^'^^ which satisfies the relation 

M^:l = (m,)^^. (10) 

If the heavy Majorana neutrino states are appreciably mixed with the active neutrino 
flavors, this equality no longer holds and LNV processes do not provide direct limits on 
Majorana neutrino mass matrix elements. 

From the non-observation of the LNV processes one can deduce the upper limits on 
the corresponding parameters {m,^) and (M^^). It must be stressed that these limits have 
physical sense only if they satisfy the following consistency conditions 

|(m^)a/3| <go, l(M^^)a/3r^ > go, (n) 

which follow from the conditions of Eq. © • 

Currently the most stringent limits of this type stem from the Oi//3/3-decay. Its amplitude, 
written in the form of Eq. ((7j), depends on the parameters {mj)ee and (M~^)ee- Assuming 
that only light or heavy exchange mechanism is in operation, the following limits have been 
derived from the experimental data 



|(m^)ee| < 0.55 eV, | (M^^)ee| ^ > 9 x lO^GeV. (12) 

Note that these limits satisfy the consistency conditions in Eq. (jlip since the characteristic 
energy scale of Oi//3/5-decay is of the order of go ~ 100 MeV. 

As we shall demonstrate, the current and near future experimental searches for (/i^,e^) 
conversion are unable to reach meaningful limits on the corresponding parameters (m^)^e 
and (M~^)^e satisfying the consistency conditions in Eq. (Illj). Moreover, the limits fol- 
lowing from the neutrino observations and cosmological data show that the sensitivities of 
(/i^, e"*") conversion experiments are too far from being able to detect neutrino contributions. 
With the lucky exception of the Oz//?/5-decay this is the fate of all the experiments searching 



for other known LNV processes (see, for instance, 42| ) 



III. EFFECTIVE NEUTRINO MASS FROM NEUTRINO OBSERVATIONS 

Here, we estimate the effective light {m^)f^e and heavy (M^^)^e neutrino effective masses 
which determine light and heavy Majorana neutrino contributions to (/i~ — e"^) conversion 



according to the general formula in Eq. (|7j). To this end we utilize the existing neutrino 
oscillation, cosmological and accelerator data, applying the methods previously used for the 
analysis of {m^)^^ relevant for Oi//?/?-decay (see, for instance, and references therein). 

Let us start with the three light neutrino scenario without heavy neutrinos. In this case 
we have 

I {mi,)^e\ = \UelU^imi + f/e2^//.2"^2 + UesU^sTJlsl , (13) 

with the unitary Pontecorvo-Maki-Nakagawa-Sakata neutrino mixing matrix U. In its stan- 
dard parametrization (e.g. [s^) it takes the form 



U 



where c, 



C12C13 



■S12C13 



'■512C23 ~ Ci2S23Si3e''^ C12C23 — Si2S23Si3e'^ 



Si3e 
■523C13 



■A 



S12S23 - Cl2C23Sl3e"^ -C12S23 - Si2C23Si3e"^ C23C13 



1 



&^ 











e 2 



(14) 



cosft 



sin 6ij. The three mixing angles vary in the range < 6ij < 7r/2. 



In addition, Majorana neutrino mixing matrix U contains three CP-violating phases: one 
Dirac 6 and two Majorana phases 0:21, 031. 

The global analysis of the solar, atmospheric, reactor and accelerator neutrino oscillation 
data gives the following values of the neutrino mixing angles 



sin^^is = 0.30 [0.23 -0.39] 
sin^ 9i3 = 0.006 [< 0.054] 
sin^ ^23 = 0.52 [0.31 - 0.72] 

and the two independent mass-squared differences ^. 



(15) 
(16) 
(17) 



Ami, = 6.9 X 10^^ eV^ [(5.4 - 9.5) x 10"^ eV^] 

-3 „-\r2 r/1 A Q '7^ 1 n-3 „-^A2l 



''sol — 

Aml^^ = 2.6 X 10"^ eV^ [(1.4 - 3.7) x 10"^ eV 



(18) 
(19) 



The values in the square brackets correspond to the 3a intervals. 

Using the above best values for the neutrino oscillation parameters we estimate the ef- 
fective light Majorana neutrino mass \{m,^)^e\ for the three standard cases of neutrino mass 
spectrum. 



^ Mass-squared difference is defined as Am^- = — 
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(1). Normal hierarchy: mi <^ m2 ^ m^. In this case Am'21 ~ Am'j^i, ~ ^"^ttm 

Therefore, one has 



mi < \ Aml^i, m2 ~ A/Am^^,, mg ~ V Am^^^. (20) 



(2). Inverted hierarchy: <^ mi < m2. Now, Am^i ~ Am^^i, Aml^ ~ — Am^^^. This 
results in the following estimate for neutrino masses 



< \/Am2j^, m2 ~ ^/Amlt^, mi ~ v^Am^j^. (21) 



Using the estimates (f^mi - ipTj) in Eq. (fT^ with the best-fit values for the neutrino oscilla- 
tion parameters from Eqs. (P3j) -()19 p . we end up with the values of the effective light neutrino 
mass for 

Normal hierarchy : \{m^)f,e\ — (0.35 - 5.3) x 10"^ eV (22) 
Inverted hierarchy : \{m,y)^e\ — (0.3 — 3.3) x 10"^ eV. (23) 

within the ranges corresponding to the variation of CP-violating phases within the intervals 
< 6 < 2tt, < ai2 < 2tt, < 0:23 < 27r. The small terms with mi in Eq. (j22j) and m^ in 
Eq. (|23p were neglected. The effect of these terms is presented in Fig. ^ which shows the 
dependence of the allowed regions of \{m^)f^e\ on the mass of the lightest neutrino mi for 
the normal and m^ for the inverted neutrino mass hierarchies. 

(3). Quasi- degenerate hierarchy: mi ~ m2 — m^. This mass spectrum can be consistent 
with neutrino oscillation data if the characteristic neutrino mass scale is sufficiently large 

mo ^ a/ Am^j^. In this case the effective light neutrino mass can be written as 

3 

lY^U^^kUekl (24) 

k=l 

In order to estimate its value one needs the values of the characteristic neutrino mass scale 
mo. It can be deduced from '^H experiments and cosmological data. Using the best fit values 
of neutrino mixing angles from Eq. ()15|) and adopting for the simplicity 6 = ai2 = 0:23 = 
we obtain 



|(^i/)^te| ^1-46 eV, mo < 2.05 eV, Troitsk experiment 
\{fnu)fj.e\ ^1-56 eV, mo < 2.2 eV, Mainz experiment 



\{m^)f,e\ <0.16 eV, mo < 0.23 eV, Cosmological data 4^| (27) 



\{fnu)fj.e\ ~0.14 eV, mo ~ 0.2 eV, Cosmological data 



4J (25) 
(26) 



rT 



m (28) 
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Note that the results of the global analysis of the cosmological data in Refs. 

provide significantly more stringent limits on the neutrino mass scale than those from the 



direct laboratory measurements of /3-decay j4J], 45|. However, at the same time the 
cosmological limits are more model dependent than the laboratory ones. 

Now, let us assume that there exist heavy neutrinos with the masses ^ go ~ ^fi, 
where go ~ f^ij. is the typical energy scale of — e"^) conversion set by the muon mass m^. 
Their contribution to this process is determined by the effective mass 

{M^').e= E ^ (29) 

fc=heavy 

Due to the lack of model independent information on mixing matrix elements U^kUek in the 
sector of heavy neutrinos it is hard to estimate this quantity. For this reason we adopt the 
conservative upper bound following from the existing LEP limit on the mass of heavy stable 
neutral lepton M^r > 39.5 GeV j4y|. Assuming the existence of only one heavy neutrino 
identified with this particle we obtain 

|(M^')^e| < (39.5 GeY)-\ (30) 

In what follows we will use the results presented in Eqs. fl2^ . (j^ . and (jHU)) for 

discussion of the expected rates of (/x~ — e+) conversion induced by the Majorana neutrino 
exchange. 



IV. NEUTRINO MEDIATED (At-,e+) CONVERSION. GENERAL FORMALISM 

The process of (/i~,e"'") conversion is very similar to the Oz//5/5-decay. Both processes 
violate lepton number by two units and, therefore, take place if and only if neutrinos are 
Majorana particles with non-zero mass. 

On the other hand, there are various important differences between (/i~,e"'") conversion 
and Oz//?/3-decay. Among them we mention the following. 

i) They have rather different available energies and different number of leptons in their 
final states. This results in a significant difference between the corresponding phase 
space integrals. 

ii) The emitted positron in (/i~,e"'") conversion has large momentum and, therefore, the 
long- wave approximation is not valid in contrast to Oi//3/3-decay. 

10 



iii) As we will show, the nuclear matrix element of (/i~,e^) conversion for light neutrino- 
exchange demonstrates a singular behavior, absent in the Oz//5/5-decay. This feature 
gives rise to the large imaginary part of the {fi~, e'^) conversion amplitude. Technically 
the singularity significantly complicates the numerical calculation of the nuclear matrix 
elements. 

iv) In the case of the e"*") conversion there is large number of nuclear final states which 
must be properly taken into account. 

Below, we analyze the amplitude of the {fi~,e~^) conversion in nuclei mediated by light 
and heavy Majorana neutrinos. The corresponding diagrams are shown in Fig. |21 We 
concentrate only on the nuclear transition connecting the ground states (g-s) of the initial 
and final nuclei, which is favored from the experimental point of view due to the minimal 
background. The characteristic signature of g.s g.s. transition is the presence of a peak 
in the spectrum at the energy 

E,+ =m^-6^- {Ef - Ei) (31) 

which allows reliable separation of signal from background. Here, m^, Sb, Ei and Ef are 
the mass of muon, the muon atomic binding energy (for ^^Ti this is Eb = 1.45 MeV), the 
energies of initial and final nuclear ground states, respectively. Latter on we neglect the 
kinetic energy of final nucleus. 

The leading order (/i~,e^) conversion matrix element, corresponding to the diagrams in 
Fig. 12 reads 



X 



(32) 

^e^A Ue_^(,e+)* ^ 2nS{E,- +E,-Ef-E,,). 



2nR L 

Here mg and are electron and proton masses, kf,+ {Ee+), k^- [E^-) are the momentum 
(energy) of outgoing positron and captured muon respectively. The conventional normaliza- 
tion factor involves the nuclear radius i? = 1.1 A^^^ fm. For the weak axial coupling constant 
gA we adopt the value gA = 1.254. In the above expression we introduced for convenience 
the following LNV parameters 

<^ = ^^, < = (M^i),,m,. (33) 
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The nuclear matrix elements in Eq. ()32|) defined as 

_^(^e+)*^ ^ + M(f4)* for t = u,N (34) 

contain the Fermi m'^^ and Gamow- Teller Mq'^ contributions. They take the follow- 
ing form for the light Majorana neutrino exchange mechanism 



+ 



(0.^1 Y.m'r+e^'^-''-^{r^)\n){n\ r+e-^^^e+-ne-ig-n |o+) ' 

g + i?e+ + En - Ei+ iSn 



nApe+)'S> _ 47ri? f dq ^ 2^ 



X 



H : 

q + Ef,+ + En — Ei + iSn 

and for the heavy Majorana neutrino exchange mechanism 

(37) 

with 

hF{q') = fW), hoT{q') = Si ■ amfliq'). (38) 
We use the conventional dipole parametrization for the nucleon form factors 0| 

/v(g ^) = (i + ^) ' ' /A(g ^) = (i + IJ) ' ' (39) 

with Av = 0.71 GeV, Aa = 1.09 GeV. In Eqs. (j2Sl)-(inZI) the factor ^{r) is the radial part 
of the bound muon IS wave function (see Appendix El) • In the denominators of Eqs. 
(jnni) we introduced the widths of intermediate nuclear states. 

In the calculations of nuclear matrix elements we adopt the following approximations. 
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i) Taking into account slow variation of muon wave function within the nucleus we apply 
the standard approximation jlQ^ 



\M 



(/^e+)$|2 



2| \y((fJ-e^)\2 



(40) 



Here ($) is the muon average probability density and 



M 



(^ie+)<I> 



$=1 



(41) 



The explicit form of ("l*)^ is given in Appendix iBl 



ii) In muon to positron conversion the typical energy of light intermediate neutrinos is 
about 100 MeV {uj \q\ > 1/ R 100 MeV) which is much larger than the typical 
excitation energies of intermediate nuclear states. Therefore, to a good approximation 
the individual energies of these states in the energy denominators of Eqs. (j35|) . (j36|) can 
be neglected or replaced by some average value < > to which the matrix elements 



are not very sensitive. Then the intermediate nuclear states can 
closure. A similar situation occurs in the case of Oi//5/3-decay 



3e summed up by 



Mini- 



Thus, in Eqs. fl35p . (jHBj) we complete the sum over the virtual intermediate nuclear 
states by closure after replacing e„ with some average values {En), e, respectively: 

\n){n\ 1 



n 

E 



q-E^-+En-E, + ien q-E^- + -Ei + ie' 

\n){n\ 1 
q + E,+ + En- Ei + ien^ q + + {En) -Ei + ie' 



(42) 
(43) 



Obviously, the validity of the closure approximation is just the question of the choice of the 
average excitation energy which will be discussed in Section [V] 

The angular part of neutrino propagators can be integrated using the relation 

-ig-(n-r™)g-ife^+.ri^^^ = 
(47r)2 5^(-l)^v/2ATTjA(fce+i?Zm)jo(gn^)jA(fce+n„/2) {Yxi^rJ ® Ya(1)k,„)}oo , (44) 



Where Ja is the spherical Bessel function, is the spherical harmonic and 

fi + fj 



ly ^ ty fy /y 



Rij I -Rij I . 



(45) 



13 



Note that in the hmit when the outgoing positron momentum \ke+\ is zero the right hand 
side of Eq. (jl^ is reduced to 47rjo(grjm). 

With the above approximations and comments we can write down the expressions for the 
nuclear matrix elements introduced in Eq. (PT|l in the form 

Here the nuclear matrix element Aii^"^ ■* is decomposed into the contributions coming from 
direct and cross Feynman diagrams in Fig El They can be written as 

^cSrf ^ = (0+1 $^r+r+47r5^(-l)^v/2ATIjA(A:e+i?/™)jA {Yxi^nJ^Yxi^RJU 

'a / 

(47) 



X-/ — ^ Z/f7 F^> r^-^-™)-^h^^'Q/)' 

Jo 1- E^- + {En) -Ei + \e \ gX J ^ 



M^:^ = (0+1 J]r+r+47r J](-l)^v/2ATTjA(A:e+i?.„.)jA {YaI^^., J ® Y,(fi«, J} 

X-/ ^ -TTTT ^ , ■ {(^l- ^m/A g ^Uc^gOJ. 

TT Jo g + + (K) - + 16 V 5-1 / 



^ Im 



— {0t\J2r+T^4nJ2i-iVy^^y+^ix{ke+Rim)}x (- 



(48) 

The Gamow- Teller and Fermi nuclear matrix elements of heavy Majorana neutrino ex- 
change mechanism take the form 

'K+rir_ 

(49) 

r) C) POO ^ ' 

x{YA(a,„)®YA(fi«,„)}oo— / io{qriMq')q'dq\Op (/ = F,GT), 

^ Jo 

with hi{q'^) defined in Eq. 

It is important to note that the value of Er = —E^^^- + {E„)—Ei is negative for the studied 
nuclear system A = 48. Therefore, the contribution of direct Feynman diagram in Fig. 2a 
with the light intermediate neutrino has the pole at g = —Er — ie, as it follows from the 
formula ()47|). As a consequence, the imaginary part of the {fi~,e~^) conversion amplitude 
for the case of the light neutrino exchange can be significant. This fact was first noticed in 



Ref. U3 and then in Refs. 



23 



2J. In Ref. 



. 3i 



it was shown that the imaginary part of the 
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amplitude dominates in the total branching ratio of the {fi , e"*") conversion in 2^A1. In the 
next section we will demonstrate that the similar conclusion is valid for (//", e^) conversion 
in ^^Ti. 

The following comment is in order. In the expressions (j35p - (j37|) for nuclear matrix ele- 
ments Ai\^'^ ■** we neglected the contributions of the higher order terms of nucleon current 
(weak-rnagnetism, induced pseudoscalar coupling). As suggested by the analogy with Ou^P- 
decay j50|, these terms should not be essential for the light neutrino exchange mechanism 
meanwhile their contribution in the case of heavy Majorana neutrino exchange might be 
significant. However, the detailed study of this effect is beyond the scope of this paper and 
will be considered elsewhere. 

Now we are ready to write down the expression for g.s. —>■ g.s. {fi^,e^) conversion 
rate. For simplicity we assume that only one mechanism is in operation and present the 
corresponding rates for light and heavy Majorana neutrino exchange mechanisms separately: 



(A'e+)|2i (A'e)|2 



TT 



where c 



■fie 



2GUim,m^)/{A7rm^R)]'gl k, 



\k 



F{Z, E) in Eq. ()50|) we take in the standard form 

2 



r^. T, (^ = z/,A^) (50) 
The relativistic Coulomb factor 



F{Z,E) 



r(27i + 1) 



(2pi?)2(^-i)|r(7i-iy)pe--^ 



(51) 



where 71 = a/1 — {aZy, a is the fine structure constant, y = aZE/p. 

To conclude this section we point out that in our analysis of {fi^, e^) conversion we limit 
ourselves by the 0^^ —>■ 0^^ transition which represents a particular contribution to the 
total rate of this process. This is the most favored channel for experimental study since its 
signal can be reliably separated from the background as we commented above. On the other 
hand in Ref. {2^ it was demonstrated that 0^^ 0+^ transition constitutes about 41% of 
the total (/i~, e^) conversion rate in ^''Al and, therefore neglecting the excited final states is 
a reasonable approximation. We expect that this conclusion holds for ^^Ti as well. 



V. NUCLEAR MATRIX ELEMENTS 



We calculate the (/x , e"*") conversion nuclear matrix elements within the proton- neutron 
renormalized Quasiparticle Random Phase Approximation (pn-RQRPA) 34 , 



35 



51 



In 
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the present study we focus on "^^Ti nucleus utilized as a stopping target in SINDRUM II 
and PRIME js^ experiments. 

Nuclear transition scheme for the studied A = 48 nuclear system is shown in Fig. El Our 
nuclear structure calculations involve the single-particle model space both for protons and 
neutrons consisting of the full — 3hu shells plus 2si/2, 0(?7/2 and Ogg/2 levels. The single 
particle energies were obtained using the Coulomb-corrected Woods-Saxon potential. The 
two-body G-matrix elements were calculated from the Bonn one-boson exchange potential 
on the basis of the Brueckner theory. Since the considered model space is finite the pair- 
ing interactions have been adjusted to fit the empirical pairing gaps j^. In addition, we 
renormalize the particle-particle and particle-hole channels of the G-matrix interaction of 
the nuclear Hamiltonian H by introducing the parameters Qpp and Qph, respectively. The 
two-nucleon correlation effect has been taken into account in a standard way by multiplying 
the operators with the square of the correlation Jastrow-like function js^. The details of 
our nuclear model can be found in Appendix O 

As we already commented in section IIVI the matrix element of the direct contribution 
(Fig. 2a) of light neutrino exchange mechanism contains an imaginary part which stems 
from the pole of the integrand in Eq. ()47|) at g = —E^ — is. Taking into account that the 
widths e of low lying nuclear states are negligible in comparison with their energies one can 
separate the imaginary and real parts of this matrix element using the well known formula 

^ V-- m5{a) (52) 



a + ie a 

valid in the limit £ — > 0. 

In Table |l] we show the nuclear matrix elements of light and heavy Majorana neutrino 
exchange mechanisms of the (/x~,e+) conversion in ^^Ti calculated for Qpp = 1.0 and Qpp = 
0.8, 1.0, 1.2. All the presented results were obtained for the particular value of energy 
difference {En) — Ei = 10 MeV. This choice is justified by weak dependence of the matrix 
elements on this parameter within the interval of its reasonable values 2 MeV < {{En)—Ei) < 
15 MeV. We verified this property by the direct numerical analysis. In Fig. ^we present the 
absolute value of the light neutrino exchange nuclear matrix element IM*^^^^-*! as a function 
of the average value (En) — Ei for Qpp = 0.8, 1.0 and 1.2. One can see that its variation 
within the studied range of {En) — E^ is about 30%. For Qpp = 0.8, 1.0 {gpp = 1.2) the 
matrix element is an increasing (decreasing) function of {En) — E^. Different behavior in 
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these two cases is related to a specific interplay between the direct Mj^^^ ^ and cross AdcTo. 
diagram terms in M^^^'^\ For gpp = 0.8, 1.0 there is a mutual cancellation of the real parts 
of these two terms so that the imaginary part of Mj(^^ which is a growing function of 
(En) — Ei, dominantes and determines the behavior of M^'^^^^ For gpp = 1.2 the situation is 
opposite. The real parts, decreasing with (En) —Ei, contribute coherently and constitute the 
dominant part of M*-^^^-* which becomes a decreasing function of (En) — Ei. We have also 
found that the nuclear matrix elements do not show an appreciable variation in the physical 
region of the parameter Qph (0.8 < Qph < 1.2). On the contrary, as seen from Table |l] they 
significantly depend on the renormalization parameter Qpp and on the two-nucleon short- 
range correlation. It is also worth noting that the large momentum /Ce+ of outgoing positron 
is the source of strong suppression of the (/i~,e"^) conversion matrix elements. In order 
to illustrate this effect we presented in Table IJ the matrix elements calculated in the limit 
|fce+| = when the suppression of this type is absent. The cross check of Table H] reveals the 
corresponding suppression factor of about ~ 10. 

An important issue of our analysis is the presence of the significant imaginary part of 
matrix element Ai'f,'^^ corresponding to the light Maiorana neutrino exchange mechanism. 

n nn 

This fact was first noticed in Ref . and then in Refs. ^31 UM ■ In the previous studies of 

n n n n 

(/i , e"*") conversion ll^, lla, the role of imaginary part was overlooked. 

In the presented detailed study we have found, that the relative contribution of the 
imaginary part to the rate of {fi~, e'^) conversion in ^^Ti is always significant but appreciably 
depends on the value of the nuclear model parameter Qpp and on the short range correlations. 
It absolutely dominates over the real part by the factor of ~ 16 for the most conventional 



case when Qpp = 1 and the short range corre^ 
of this choice see, for instance, Ref. [3, 



ations are taken into account (for the motivation 
IJ]). This conclusion is consistent with the result 
of Ref. studying (/i~,e^) conversion in ^^Al within shell-model approach where it was 
found that the imaginary part for the light neutrino exchange dominates over the real one by 
the factor of about 20. However it is notable that the relative contribution of the imaginary 
part is model dependent and can vary from one nucleus to another. In this situation the 
role of the imaginary part in (/i~,e+) conversion requires further study for other nuclear 
systems. 

From the view point of nuclear structure theory it is instructive to compare the values 
of (/x~, e'^) conversion nuclear matrix elements with the corresponding values of Oz//?/3-decay 
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matrix elements of A = 48 nuclear system. For Oi//3/3-decay this system is represented by 
^^Ca with the matrix elements 

IM'i^^^l = 0.82, |A^5^'^| = 24.2 (53) 

derived within the pn-RQRPA approach in Ref. 14 • As seen, the matrix elements of the 
(/i~,e"'') conversion (fS3j) are strongly suppressed in comparison with those of Oz//?/5-decay 
by the factors of about 17 and 5 for the light and heavy Majorana neutrino exchange 
mechanisms respectively. As we commented above the explanation of this difference between 
the two processes mostly resides in the large momentum of outgoing positron produced in 
{fi~,e~^) conversion. 



VI. (/u ,e+) CONVERSION AND EFFECTIVE NEUTRINO MASSES 

Now, let us discuss the possible issues of (/i^,e+) conversion experiments for neutrino 
physics. 

From Eq. ()50p we obtain the (/i~,e^) conversion branching ratios in ^^Ti for the light 
and heavy Majorana neutrino exchange mechanisms: 

-pipe) 

I^r ^ = ^ = 2.6 X lQ-''\Mt'\''\v\''\'' = u^N). (54) 

Here we use the known experimental value F*^^'^^ = 2.60 x 10^ s~^ [55] of ordinary muon 
capture rate in ^^Ti. For the further discussion we choose the following sample values of 
nuclear matrix elements of ^^Ti from Table IJ 

= 0.025, |A<S^"^^| = 5.2 (55) 

corresponding to Qpp =1.0 with the presence of the two-nucleon short range correlations. 
Substituting these numerical values of nuclear matrix elements to Eq. ()54|) we obtain 

Ri^'^^^ = 1.6 10-^5 X (56) 

ml 

R^^'"^ = 7.0 10-'' X \{M^'),,\'ml. (57) 

From the existing experimental upper bound in Eq. Q one obtains the following limits for 
the effective masses of light and heavy Majorana neutrinos 

|(m)^e| < 1.3 X lO'' MeV |(M^^)/,er^ > 3.3 x 10"^ MeV (58) 

18 



Obviously, these limits have no physical sense since they do not satisfy the consistency 
condition in Eq. (fTT|) with the characteristic energy scale go ~ = 105MeV of 
(/i~,e"'") conversion. Meaningful limits on the parameters {m)^^, {M^^)^e, which may have 
some impact on neutrino physics, could be reached if the {fi~,e~^) conversion experiments 
would improve their sensitivities by at least 10 orders of magnitude. Clearly, such a tremen- 
dous improvement is unrealistic for the near future experiments. 

On the other hand we can estimate the expected branching ratios of {fi~, e+) conversion 
induced by the light and heavy Majorana neutrino exchange using the estimates of {m)fj_e, 
(Mj^^)^e niade in section ITTTl from the present neutrino data. Substituting the values of these 
parameters in Eqs. (fHH|) - (f^ we obtain the following results. 
Light Majorana neutrino exchange contribution: 

i) Normal neutrino mass hierarchy, | < m >^j,e \ — (0.35 — 5.3) x 10^^ eV 

~ (0.008 - 1.7) X 10"^^ (59) 

ii) Inverted neutrino mass hierarchy, | < m >^e | — (0.3 — 3.3) x 10~^ eV 

Ri^'''^ ~ (0.05 - 6.7) X 10"'°. (60) 

iii) Quasidegenerate mass hierarchy 

Ri^""*^ <1.3 X 10"^^ (m^) < 1.46 eV, Troitsk experiment Q 



Ri^''*^ <1.5 X 10~^^ (m^) < 1.56 eV, Mainz experiment 



(61) 
(62) 

Ri^""^^ <1.6 X 10"^^ (m^) < 0.16 eV, Cosmological data 46 1 (63) 
Rl^^''^ -1.3 X 10~3^ (m^) ~ 0.14 eV, Cosmological data ^ (64) 

Let us remind that the cosmological data based limits (j63|) and (I64p . albeit more stringent, 
are more model dependent than the laboratory ones and ()62|) . 
Heavy Majorana neutrino contribution: 

R^i^"^^ < 3.8 X 10-24 (65) 

All the values of {^~,e~^) conversion branching ratio in Eq. (j59|) -(|65| ) are hopelessly low 
for being detected even in a distant future. Thus, searching for {fi~,e'^) conversion can- 
not have any direct impact on neutrino physics. On the other hand any observation of 
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(/i~,e"'") conversion at branching ratios above the hmits in Eq. ()59p -()65 p would be unam- 
biguous signal of new physics beyond the simplest extension of SM with massive Majorana 
neutrinos and would imply the presence of new interactions. 

This conclusion are in a sharp contrast with Oz//5/5-decay experiments which already 
provide an important information on neutrino properties and are expected to detect neutrino 
contribution in the near future. This is due to their unique sensitivities to Oz//?/3-decay 
signal. In order to give an impression to which extent Oi//3/3-decay experiments overcome in 
sensitivities the experiments searching for e"*") conversion let us compare, as an example, 
the rates of e"^) conversion in ^^Ti and Oz//3/3-decay of ^^Ca. To this end it is sufficient 
to consider only light Majorana neutrino exchange contributions in both cases. For the rate 
of Oz//3/3-decay we have the well known formula 

2 



r(^^)=ln2Goi 
where Gqi = 8.031 x 10"^^ year"^ ^ and 



\M 



(ee)|2 



A;=light 



(66) 



(67) 



Using the value of Oi//3/3-decay nuclear matrix element A^lf^'' from Eq. we estimate the 
ratio of {fi~, e+) conversion to Oz//5/5-decay rates: 



(Me)|2 



(ee)|2 



{rriu) 



fie 



{my). 



351 



{my) 



fie 



{my). 



(68) 



The the (/i~,e^) conversion receives a significant enhancement mostly due to the larger 
available energy of this process. Thus, for {my)^e ~ {my)ee the e"*") conversion rate rif"^ •* 
is by more than 2 orders of magnitude larger than the rate F^f^'' of Oi//5/?-decay. Nevertheless 
the experimental prospects for searching for Oz//5/5-decay are incomparably better than those 
for {fi~,e~^) conversion. This is mainly because the number of potentially Oi//5/3-decaying 
nuclei monitored in OujSlS experiments is by many orders of magnitude larger than the 
number of mesoatoms created by muon beams in the muon- conversion experiments. 



VII. SUMMARY AND OUTLOOK 



In summary, the light and heavy Majorana neutrino exchange mechanisms of , e^) 
conversion have been studied. Special emphasis was made on the nuclear structure aspects of 
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this process. We have performed the first reahstic calculations of the corresponding nuclear 
matrix elements for ^^Ti nucleus used as a stopping target in the current js^ and the 



forthcoming [33| (/i~,e"'") conversion experiments. Our analysis is based on the pn-RQRPA 
approach and limited to the case of 0^^ —>■ Oj"^ transition channel, which is most relevant 
for experimental searches for e^) conversion. The effects of the ground state and two- 
nucleon short-range correlations have been properly taken into account. We pointed out that 
their inclusion results in the significant reduction of (u~, e"*") conversion matrix elements. 

Ou. detailed analyse con«.,ned the con.ect.e of Refs. Q Q o„ the i.po.anee of the 
imaginary part of the nuclear matrix elements for the case of the light Majorana neutrino 
exchange mechanism of e"*") conversion. The similar result was recently obtained in Ref. 



24| for (/i , e"^) conversion in ^^Al. 



We also derived the limits on the effective masses of light {m)^e and heavy (M^^)^e Majo- 
rana neutrinos from the neutrino oscillations, tritium beta decay, accelerator and cosmolog- 
ical data. Using these limits we estimated the expected rates of e~^) conversion induced 
by Majorana neutrino exchange. Their values were found to be so small that even within 
a quite distant future the {^~,e~^) conversion experiments will hardly be able to detect the 
neutrino contribution and, thus, to have a direct impact on neutrino physics. On the other 
hand the eventual observation of (/x~,e+) conversion at larger rates would be unambiguous 
signal of new physics beyond the standard model implying new non-standard interactions. 
Moreover, this observation, independently of the (/i~,e"'") conversion rate, would definitely 
prove that neutrinos are Majorana particles as follows from the "black box" type theorem [l| 
establishing the fundamental relation between LNV processes and Majorana nature of neu- 
trinos. In view of this it remains actual to study possible scenarios of new physics consistent 
with the values of {fi~ , e~^) conversion rates within the reach of the present and near future 
experiments. 
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APPENDIX A: BOUND MUON WAVE-FUNCTION 



The bound muon wave function (IS-state) is given by the expression 



^{x) = C&(r)p-i^M-"° ^ (Al) 



where the radial $(r) and the spinorial parts have the forms 



73/2 

$fr) = — e"^''/"'^ (A2) 



and 



u:. = v/2l^ I I , (A3) 

with = 47r/(m^e^) [a^/a^ ~ rrie/rn^ 5 x 10^'^), is reduced mass of muon atom, Z 
is nuclear charge. 



APPENDIX B: MUON AVERAGE PROBABILITY DENSITY OVER NUCLEUS 

Muon average probability density over nucleus is defined as 

J p[x)d-^x 

where p(a;) is the nuclear charge density. To a good approximation it can be written in the 
following compact form 

= 1^, (B2) 



Here the effective charge for Z = 22 nuclear system is Zefr is Z^s = 17.5 



0. 



APPENDIX C: NUCLEAR MODEL 



Here we shortly outline our approach to the nuclear structure calculations. 

We introduce particle (quasiparticle) creation operators as cl^^ («tmr) ^'^^ t = p,n. The 
indices p = {np,lp,jp) and n = {nn,ln,jn) denote proton and neutron quantum numbers in 
a particular shell. Transformation from the particle to quasiparticle basis is realized by the 
Bogolyubov transformation 
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where the tilde denotes time reversal, drm^ = (— l)•'^~™^aT-mT• 
Occupation amplitudes Ur, Vr and quasiparticle energies Er are obtained by solving BCS 
equation jsTj 

-Sr + XrJ yVrJ \Vr J 

where Er is the energy of single particle state derived from the Wood-Saxon potential. The 
pairing potential takes the form 

A. = (2jV + 1)-'/' J2^2ja + lY^'Giaa, tt; J = 0)m,z;,. (C3) 

a 

Here G{aa,TT; J) is particle-particle matrix element defined e.g. in Ref. 3|. The value of 
Lagrange multiplier A is fixed by the particle number in non-correlated BCS vacuum 

(iV,) = 5^(2j. + l)vl (C4) 

T 

After the diagonalization, BCS Eq. (j('2jl reads 



Er = ViSr - Xr? + A^, = ^ (1 _ £l_A), ^1 = 1- V^. (C5) 

This system of equations can be solved by the iteration of the parameter with the con- 
dition = (A^^). 

The nuclear Hamiltonian in quasiparticle representation takes after the BCS transforma- 
tion the form 

H = ^ ErO^^^armr + H22 + H40 + HqA + -f^31 + -^13, (C6) 
rrriT- 

where Hij is the normally ordered part of residual interaction with i creation and j annihi- 
lation operators. 

Within pn-RQRPA, the m-th nuclear excited state |m, JM) with the angular momentum 
J and its projection M is obtained from the RPA vacuum lOj^p^) 

\m,JM)=Qf^.\Qi^^), (C7) 

where RPA vacuum is defined by the condition 

QjM- 1 Orpa) = 0- (C8) 
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and phonon operator Q'J'm^ is defined as 

Qjil^ = ^^nJ-^)^\'on.JM\ ~ ^(^n,J^)\pn,.JM)-, (C9) 

pn 

where Aj^^ (A(p„ jtt)) is two-particle creation (annihilation) operator which couples quasi- 
particles to the angular momentum J with the projection M: 

mp,m„ 

A{pn, JM) = (-l)-^-^A(pn, JM) = (-l)-^"^ J] Cj-J^.v^n^^Sn.. (dl) 

Here C/^ , ™ are Clebsh-Gordan coefficients. 

The commutator [A, A^] is replaced within pn-RQRPA by its mean value in the QRPA 
vacuum 

[AA^] - (0+pA|[A(pn, JM),A(pV, JM)]|0+pa) 

= Spp,6nn' I 1 - -^{OtpAlH^phoKpA) " ^ (OrPA I [4an]oO |0^pa) > (C12) 
Jp Jn J 

where jp = ^/^jp + 1 and 



Within the quasiboson approximation, RPA vacuum lO^p^) in Eq. ()C12|) is replaced by non- 
correlated BCS vacuum lOgcg) (i.e. Dp^^j^ = 1). Quasiboson approximation violates Pauli 
exclusion principle. 

From the Schrodinger equation 

[H, Q^M'r] |Orpa) = ^J'^'^jM'^ I Orpa)' (C14) 
with the excitation energy flj-n, we obtain RQRPA equation, 
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Here matrices A, B have the form 

•^pn,p'n' = i^P + En)5pp'5nn' - 2[G{pn,p'n' \ J){UpUnUpmn> + Vp'VnVp'Vn' 
+ F{pn,pn'; J){UpVnUpfVn' + VpUnVp'Un')]Dl^^j^Dl!l j^, 



^pn,p'n' = i^P + En)2[G{pn,p'n'; J){UpUnVp'Vn' + VptVnUpiUn') 

- F{pn,pn'; J){UpVnVp>Un' + VpUnUp>Vn')]D]l^^j^D^Jl,^j^ 

TfL TTh 

and amphtudes X, J^^, X, are 



(C16) 



(C17) 



^ipn,J-^) — ^pn,J^^(pn,J^)i (pn,J^) — ^ pn,J^ {pn,J^)i V^^^) 

where F{pn,p'n'] J) is the particle-hole interaction matrix element. From the mapping 
procedure ()C12|) we obtain for the coefficients Dpn^j the system of nonlinear equations 



1 ^ 1^ -^P"'..^'"l^{pn',J-)P ~ ^ y2 ^P'^^J"'\^(p'n,J^)\^- (C19) 

p n'J'm Jn p'J'm 



The amplitudes j^^, l^(p„ j^r) and the excitation energies ^'J'^ are obtained by iterating 
of the coupled equations ()C19|) a ()C15|) . 

The {fi'jC'^) conversion nuclear matrix elements within pn-RQRPA are transformed to 
the sum of the two-particle matrix elements 



Jn jp ^ 



(C20) 

X (p(l),p'(2);:r|/(ri2)r+r+Of/7(ri2)|^(l),^'(2);:r) 
X (0;||[cJ,,c„,]j||J-m;)(J-m/|J-mi)(J-m,||[c;,c„]j||0+). 

Here {• ■ ■ } is the Wigner 6j symbol, (9^2^''^ is space- and spin-dependent part of the matrix 
element. The single particle densities are defined as 



V2J + 1 



— l^^p ^(pn,J^)+^p ^ {pn,J^))\J ^pn,J-^i (^L/ZZJ 
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where the indices (i) and (/) indicate that the excitations are defined with the respect to 
the ground state of the initial and final nucleus respectively. When these states are not the 
same, the overlap factor 

pn 

must be introduced 0]. Repulsion between the nucleons at short distances is described by 
the short-range correlation factor f{r 12) of the form 

/(n2) = l-e-°'--(l-6r?2), (C24) 



where a = 1.1 fm^ a b = 0.68 fm^ 



Particle-particle and particle-hole channels of the 



nuclear Hamiltonian are renormalized by the parameters Qpp and Qph- 



F{pn,p'n; J) gphF{pn,p'n'] J), (C25) 
G{pn,p'n; J) gppG{pn,pn'] J). (C26) 
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FIGURES 
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FIG. 1: Allowed regions of the effective Major ana neutrino mass |(m)^e| for normal (left panel) and 
inverted (right panel) hierarchy vs. the mass of lightest neutrino state: mi and ma, respectively. 
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FIG. 2: Direct (a) and cross (b) Feynman diagrams of (// ,e'^) conversion in nuclei mediated by 
Major ana neutrinos. 
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FIG. 3: Transition scheme for the A = AS nuclear system. 
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FIG. 4: The nuclear matrix elements of the light Majorana neutrino exchange mechanisms of the 
(/Lt~, e"*") conversion in ^Ti as a function of the average value of energy difference < En> —Ei. 
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TABLES 
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TABLE I: Nuclear matrix elements of the light and heavy Majorana neutrino exchange mechanisms 
of {fi^,e^) conversion in '^^Ti [see Eqs. ()46() - ()49|) ]. The calculations have been performed within 
pn-RQRPA without and with the inclusion of two-nucleon short-range correlations (s.r.c). 
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